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Abstract 

Contact-drop dispensing is central to many small-scale applications, 
such as direct-scanning probe lithography and micromachined fountain- 
pen techniques. Accurate and controllable dispensing required for nanometer- 
resolved surface patterning hinges on the stability and breakup of liquid 
bridges. Here, we analytically study the stability of catenoids pinned at 
one contact line with the other free to move on a substrate subject to ax- 
isymmetric and non-axisymmetric perturbations. We apply a variational 
formulation to derive the corresponding stability criteria. The maximal 
stability region and stability region are represented in the favourable and 
canonical phase diagrams, providing a complete description of catenoid 
equilibrium and stability. All catenoids are stable with respect to non- 
axisymmetric perturbations. For a fixed contact angle, there exists a 
critical volume below which catenoids are unstable to axisymmetric per¬ 
turbations. Equilibrium solution multiplicity is discussed in detail, and we 
elucidate how geometrical symmetry is reflected in the maximal stability 
and stability regions. 


1 Introduction 


The celebrated treatise of Plateau [l| was a key study in the nineteen cen¬ 
tury, addressing liquid-bridge stability and breakup. Early investigations were 


, crystal growth in 
Recently, interests 


motivated by applications such as liquid-jet breakup 
microgravity [3] , oil recovery Q , and paper wet st reng th 
have grown into areas such as elastocapillarity [THlOj, contact-drop dispensing 
111 with applications to scanning-probe lithography [l^ and micromachined 
fountain-pen techniques [l^. Molecular-resolution surface patterning provides 
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new opportunities for advanced tissue engineering Il4. DNA self-assembled 
nanoconstructs 0 , and highly sensitive protein chips [l^. 

Static stability analysis of liquid bridges can be traced to the nineteenth 
century Bin. Howe’s variational formulation extended Plateau’s primitive 
theory for cylindrical interfaces to unbounded axisymmetric capillary surfaces 
subject to a constant-volume. His criteria (sufficient conditions for the weak 
extrema of a functional) guarantee a surface to have the minimum e nerg y among 
all the neighbouring surfaces of revolution. Gillette and Dyson [ISj applied 
Howe’s method to predict the stability limit of bounded axisymmetric liquid 
bridges with respect to axisymmetric perturbations. These criteria were later 
generalized for arbitrary interfaces with arbitrary perturbations [l9l |. 

Catenoids are doubly-connected surfaces of revolution with zero mean cur¬ 
vature. They are special cases of constant-mean-curvature axisymmetric sur¬ 
faces, and are important to stability studies on weightless liquid bridges for 
two reasons: (1) Stability criteria can be obtained analytically for catenoids, 
which helps guide numerical algorithms for general liquid brides in the small 
pressure (mean curvature) limit, and (2) the curve corresponding to catenoidal 
interfaces in the volume-slenderness phase diagram defines a boundary between 
regions of positive and negative capillary pressure [l9l |. This is important for 
mechanical systems wit h dy namics that are driven by capillary pressure {e.g., 
elasto-capillary systems [23). Previous studies have considered catenoids bridg¬ 
ing two circular discs of the same radius 121|. c atenoids between two parallel 
plates with both contact lines free to move [221424 1. and catenoids between a 
plat and sphere 0. However, these results are not applicable to contact-drop 
dispensing applications where the liquid forms a bridge with a free contact line 
at one end. 

Recent studies on contact-drop dispensing have shown that the deposited 
drop size can be adjusted by the needle retraction speed, needle-tip size, sur¬ 


face characteristics, and dispensing control parameters [111.I26j. However, these 
studies do not distinguish the effect of static parameters from dynamic ones. A 
static stability analysis of liquid bridges with a free contact line furnishes the 
maximum-height stability limit, which reasonably approximates the pinch-off 
height at small capillary numbers Ca <C 1 (quasi-static limit) 0,0 ■ Critical 
perturbations estimate the dispensed drop volume and show how the bridge dy¬ 
namically evolves. Here, we only focus on the catenoid as an important special 
case since equilibrium solution multiplicity and stability criteria can be deter¬ 
mined analytically. 

In this paper, we analytically study the static stability of catenoids pinned 
at one contact line with the other free to move on a flat substrate with respect 
to constant-volume perturbations. This furnishes a two-dimensional phase di¬ 
agram in which the stability region is represented with respect to the catenoid 
volume and slenderness. The effect of the catenoid geometrical symmetry on 
the stability region boundaries is discussed. We also present a phase diagram 
with respect to canonical variables, which facilitates the representation of sym¬ 
metry in the stability region, maximal stability region, and multiple equilibrium 
solutions subject to various constraints. Myshkis et al. 0 described how free 
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Figure 1: Catenoidal liquid bridge: schematic (left) and coordinate system with 
meridian curve parametrization (right). 


contact lines are generally treated in their variational method. However, the 
stability criteria were not presented for liquid bridges with free contact lines. 
Therefore, we first present an exposition of Myshkis’s variational formulation 
[l^ , and then derive the stability criteria in section [ 2 ] for axisymmetric liquid 
bridges with a free contact line. Equilibrium solution multiplicity is discussed 
in section o and the maximal stability and stability regions are determined 
for cylinders and catenoids in sections [3. 2. II and 13. 2. 21 respectively. The results 
are summarized in sections |4l 


2 Theory 

We consider a liquid of volume v bridging a circular disk with radius Rq and a 
large plate. The disc and plate are separated by a distance h, as shown in Fig.[TJ 
The region occupied by the liquid bridge is denoted fli, and that occupied by 
the surrounding fluid (of a different phase) is denoted fig. The bridge is pinned 
to the disc and free to slide horizontally on the plate. We restrict our analysis 
to catenoidal liquid bridges, which implies that the regions g and I have the 
same density and pressure, and the interface Tgi is a surface of revolution. The 
formulation is presented as the limit of axisymmetric weightless liquid bridges 
with mean curvature approaching zero. The origin of the coordinate system 
is placed on the plane passing through the catenoid neck such that the z-axis 
is the symmetry axis. The meridian curve is parametrized with respect to its 
arclength s such that s = 0 at z = 0. An equilibrium surface is specified by 

se[so,si], (1) 

extremizing the potential energy 



U = 'yslTsl + Iglfgl + Isgfsg, (2) 

where 7 ^- is the surface tension between the phases i and j, and F^- is the 
interfacial surface area. Following Myshkis et al. [l^, this leads to the Young- 
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Laplace equation 


r r" = -z'{q-z'/r) 
\ z" = r'{q-z'jr) 

for axisymmetric equilibrium surfaces and 


('= d/ds) 


7gi cos Qc. = Isg - IsU and 


cos 0c = n • Hp, 


(3) 

(4) 


where q = —2km and 9c is the contact angle. Here, km is the mean curvature, 
which is zero for catenoids. Equation Q shows that the contact angle is a 
thermodynamic property of the three-phase {g, I, and s) contact line, which is 
a constant for a specific substrate (plate) and the fluids occupying fig and fli. 
Note that the dihedral angle 6d can vary independently with the bridge volume 
to extremize the potential energy. Introducing the following lengths, which are 
scaled with the neck radius. 


= r/ro, z = zlro, s = s/ro, 


(5) 


the cylindrical volume V = v/{'KRQh) and slenderness A = are two di¬ 

mensionless parameters with which to present the phase diagram. 

Following the method of Myshkis et al. [l^, the interface stability is de¬ 
termined by the sign of the second variation. Using the Ritz method [29j, the 
second variation is associated with the eigenvalues of the corresponding Strum- 
Liouville problem. Stability studies are generally concerned with determining 
stability regions in the phase diagram. Stability-region boundaries, identified 
by 5'^U = 0, correspond to critical states, separating stable equilibrium surfaces 
from unstable ones. Hence, we seek the conditions where Aq or Ai = 0, resulting 
in 

-I- M = 0 
V^o(so) = 0, 

//o' = 0 

for axisymmetric perturbations and 


‘/5o(si) +XV?o(si) = 0 


(.if — l/f'^)Lpi = 0 

</5i(so)=0, (/ji(si)-bx(/Ji(si) = 0 


for non-axisymmetric perturbations, where 

ku cos 6c - kpi 


X = 


sin 6 c 


at £, 


( 6 ) 


(7) 


( 8 ) 


.if = 


ds^ 


U_d 
f ds 


q- 


(9) 


with q = gro, x = ’"oXi 8'iid ku, kpi the first principal curvatures of the interface 
and plate at the contact line £. The solutions of Eqs. m and o can be written 


(fois) = C'iwi(s) -I- C2W2{s) + 


( 10 ) 
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V3l(s) = C4W4(s) + (^ 5 ^ 5 ( 5 ) (11) 

for axisymmetric and non-axisymmetric perturbations, respectively 0. These 
satisfy the following differential equations and their initial conditions 

Jfwi = 0 , wi( 0 ) = 0 , u;'i(0) = l, (12) 

^u;2 = 0, ^2(0) = !, u;2(0)=0, (13) 

^W3 + l = 0, W3(0) = -1/4, u;^(0)=0, (14) 

(Jf-l/r^)w4=0, W4(0) = 0, 104(0) = 1, (15) 

(Jf - l/f^)w5 = 0, u;5(0) = l, u;5(0) = 0, (16) 


where wi, W 4 are odd and W 2 , W 3 , W 3 are even functions. Note that the initial 
conditions in Eq. (HI can be arbitrarily chosen because they do not affect 
the conditions describing the critical states of equilibrium surfaces (Eqs. (1171) 
and (EOI))- The homogeneous solution of Eq. m is obtained from a linear 
combination of wi and W 2 - Erom Eq. m, the homogeneous part of W 3 makes 
no independent contribution to the general solution of ipo ■ Therefore, the initial 
conditions for W 3 are chosen such that the general solution for W 3 comprises 
only the particular part. 

The critical state of an equilibrium surface is identified by the existence of 
a non-trivial solution for ipQ or ipi. These existence conditions can be obtained 
from Eqs. m and m as 


Wl(so) W2{so) W3{so) 

w[(si) w'^isi) w'^isi) 

fs^rwids fl"fw2ds Js^"fw3ds 

Wl(so) W 2 (so) W3(so) 

Wi(si) W 2 (si) W3(si) 

Iso fso fw3ds 


(17) 


W4(sq) 

W 5 (so) 

^4(54) 

^5(^1) 

W 4 (so) 

W 5 (so) 

W 4 (si) 

W 5 (si) 


(18) 


Here, and are the critical x corresponding to axisymmetric and non- 
axisymmetric perturbations, respectively. Note that x identifies 

a critical state. Eor a fixed Fgi , the minimum eigenvalue of the Sturm-Liouville 
problem is monotonically increasing with y. Hence, Ai > 0 for y > y* (i = 0,1). 
It follows that an equilibrium surface is unstable with respect to axisymmetric 
(non-axisymmetric) perturbations if y^ < y° (y^ > y°) when y <max{y°,y^}. 
Moreover, fixed contact lines can be represented as the limiting case of free 
contact lines when y —>• 00 . Therefore, for a fixed Tg/, A —>■ —00 as y —>■ —00 
and A —>■ as y —>• 00 ; here, A is the smallest eigenvalue of the Sturm-Liouville 
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problem, and v is the smallest eigenvalue of a similar problem with tpi = Q dX £. 
Hence, 

\<v, (19) 

implying that the stability region of capillary surfaces with free contact lines is a 
subset of the corresponding stability region for the same capillary surfaces with 
fixed contact lineo The latter is termed the maximal stability region (MSR), 
a concept introduced by Slobozhanin and Tyuptsov [s^ . The critical states are 
determined by 




Wliso) W2{so) 

Wl{si) W2{h) 

fwids fw 2 ds 


wsih) 

waisi) 


( 20 ) 




W4{So) W^iSo) 
W4{si) 


( 21 ) 


For a fixed sqj the first si along the meridian curve at which = 0 {D^ = 0) 
corresponds to a critical state of the MSR with respect to axisymmetric (non- 
axisymmetric) perturbations. Note that the MSR coincides with the stability 
region for capillary surfaces with only pinned contact lines. Moreover, deter¬ 
mining the MSR for capillary surfaces with free contact lines prior to testing 
the stability criteria given by Eqs. CZl) and (115)) is necessary, since x = and 
X = X^ generally have more than one solution. Therefore, x >max{x'^, X^} indi¬ 
cates stability only for surfaces belonging to the MSR. All equilibrium surfaces 
outside the MSR are unstable. A summary of Myshkis’s method is given by 
Bostwick and Steen (sij. 


3 Results and discussion 


3.1 Equilibrium solution 

Solving Eq. ([3]) for <7 = 0 furnishes the equilibrium meridian curve 


with 


r r(s) = Vs'^ + 1 
\ z{s) = — ln(s -I- + 1) 


A = 


x/^o + 1 


: In 


Si + x/sj -|-1\ 

So + \/sq + 1 j 


( 22 ) 


(23) 


_ Si V^s j -I- 1 - sq -|-1-|-At/sq-|-1 
2A(s2 + 1)3/2 ’ 

^What ‘the same’ means here depends on how a capillary surface is specified. For example, 
as will be discussed in section ED a catenoid, such as the one shown in Fig. ^ is uniquely 
specified by sq and si. Hence, catenoids with free contact lines are being compared to those 
with the same sq and si, but fixed at si. 
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6ld = tan ^(1/so), 


(25) 


0c = tan ^(—1/si). 


(26) 


Equations (I251) - (E51) furnish four constraints on A, V, 9c, 9d, so, and si, leaving 
two degrees of freedom. Fixing two variables, the other four and the catenoid 
geometry, as shown in Fig. [U are completely specified. Therefore, among these 
six, one can select two variables to represent stability and equilibrium data 
as two-dimensional phase diagrams. Furthermore, Eqs. (1^ and (1^ provide 
a one-to-one correspondence between (so)Si) and (6c, Od). Thus, they can be 
interchanged without affecting phase diagram characteristics, albeit (9c, 9d) are 
preferred to (sqjSi) as they vary in a finite range. In this paper, the catenoid 
equilibrium solution and stability region are presented with respect to A, V, 9c, 
and 9d subject to various constraints. 

The cylindrical volume and slenderness are the two favourable quantities 
with which stability regions have been represented in the literature [l^, 
because they can be readily measured experimentally. We refer to (A, V) as 
‘favourable parameters’ and the respective phase diagram as the ‘favourable 
phase diagram’. Moreover, Eqs. (I^ - (I^ are single-valued functions of Sq and 
si. Consequently, the left-hand-side parameters characterizing a catenoid de¬ 
scribed in Fig. [T]are uniquely specified with respect to sq and si. Thus, (sq, si) 
are more convenient for representing the MSR and stability region. We refer to 
(so,si) or (9c,9d), as ‘canonical parameters’ and the respective phase diagram 
as the ‘canonical phase diagram’. 

Identifying the existence-region boundary is necessary to properly represent 
equilibrium solutions with respect to the favourable and canonical parameters. 
This is straightforward for the canonical parameters because 9d € (tt — 9c, tt] for 
9c G [0,7r]. However, determining the existence-region boundary with respect 
to the favourable parameters is non-trivial. One may naturally suppose that 
F(so,Si,A) has a minimum for a given A. Therefore, the cylindrical volume 
given by Eq. (I2T)) is to be minimized subject to a constant slenderness A, leading 
to 


\9so)s,,a 


+ 9 



= 0 , 


(27) 



So. A 


+ 1 ? 



= 0 , 


(28) 


A(soj Si) = const.. 


(29) 


where'd is a Lagrange multiplier. Moreover, Eqs. (l27)) - (l29l) define the existence- 
region boundary in the K-V space, also furnishing a lower bound on V. 

Figure [2] shows equilibrium isocontours and existence region with respect to 
the favourable parameters. The Om curve (dashed line) is the existence-region 
boundary that corresponds to (A, V) satisfying Eqs. (ETll - d^ . No catenoid can 
be found with volume and slenderness below this curve. 
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Figure 2: Equilibrium solution and existence region (above dashed line) with 
respect to the favourable parameters; (a) constant-dihedral angle isocontours 
(labels denote 6 d in degrees), (b) constant-contact angle isocontours (labels 
denote 9c in degrees), and (c) constant-dihedral and contact angle isocontours. 


Iso-0d curves are plotted in Fig. ura), facilitating the representation of so¬ 
lution multiplicity with respect to A, V, and 6d- Isocontours can be viewed 
as the level curves of a multivalued function 9d = 9d{A.,V). For example, the 
level curves with 9d = 110° and 9d = 170° intersect at A, implying that two 
equilibrium solutions exist for the corresponding (A, E). The isocontour with 
9d = 7r/2 is a special case that asymptotes to Om. It can be proved that 


lim E(so,si)|(So,Ji)£z„ = lim E(so, si)|(so.si)GOm = lim ^ = oo, 

A—>00 A->oo A—>00 ZA 

(30) 

where Zu is the isocontour with 9d = 7r/2. Figure [2])b) can be similarly inter¬ 
preted. Here, isocontours can be viewed as the level curves of a multivalued 
function 9c = 9c{A,V). Isocontours with 9c = 150° and 9c = 170° intersect at 
B, implying that two equilibrium solutions exist for the corresponding (A, E). 
For a given 9c, there is a slenderness above which there are no catenoids. Two 
equilibrium solutions exist for smaller slendernesses. lso-9d and iso-0c curves 
are plotted together in Fig. [21(c). This figure demonstrates that isocontour in¬ 
tersections do not always correspond to equilibrium solutions. For example, the 
isocontours with 9d = 170° and 9c = 130° intersect at A, C, and D correspond¬ 
ing to three different (A, H). However, D is the only intersection representing 
an equilibrium at 9d = 170° and 9c = 130°. Note that these isocontours are 
projections of the corresponding level curves onto the (A, V) plane. Here, A and 
C lie at the intersections of the projections and do not represent equilibrium 
solutions. Hence, the space (A, V) is unsuitable for representing the stability 
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of catenoids because it cannot be partitioned into mutually exclusive regions of 
stability and instability. 

Figure [3] shows equilibrium isocontours and the existence region with respect 
to the canonical parameters. As previously discussed, equilibrium solutions for 
catenoids cannot be conveniently represented with respect to the favourable pa¬ 
rameters because isocontour intersections are not always associated with equi¬ 
librium solutions. Moreover, two solutions may correspond to the same point 
(A, V) in the existence region. In contrast, the canonical parameters furnish a 
one-to-one correspondence between points in the existence region and equilib¬ 
rium solutions. Here, the existence region is the upper triangle indicated by 
Z'Z"U. Note that the existence-region boundary Z'Z" does not correspond to 
Om in Fig. [2] Iso-A and iso-H curves are plotted in Figs. [3l[a) and (b). Iso¬ 
contours are the level curves of the single-valued functions A = A(0c, Od) and 
V = V{9c,dd) given by Eqs. (1331) and (1331) . Iso-A and iso-H curves are overlaid 
in Fig. He). Here, unlike the favourable parameters, every isocontour intersec¬ 
tion uniquely represents an equilibrium solution. Interesting to note are the 
two equilibrium solutions corresponding to the point A in Fig.H which are de¬ 
noted by A and A', and are distinctly represented with respect to the canonical 
parameters. 


3.2 Stability 


Several factors affect the equilibrium state and stability of capillary surfaces, 
including fluid inertia, external fields {e.g., gravitational and centrifugal forces), 
and boundary conditions at contact lines. The latter differentiates contact- 
drop dispensing applications from classical liquid bridge problems where the 
equilibrium surface is pinned to two coaxial discs. The contact-line condition 
can be easily accounted for in the equilibrium solution by integration constants 
of the integral curve obtained from Eq. ([S])- The influence on stability, however, 
is not straightforward. It affects the eigenvalues of the Sturm-Liouville problem 
through the boundary condition at i. This plays a far more significant role 
in the stability of capillary surfaces. For example, the notion of wavenumber 
introduced for classifying equilibrium solution branches and characterizing the 
bifurcation of liquid bridges 3^ 371 is directly related to the conditions at the 
contact lines. Before proceeding to catenoids, we elucidate the contact-line 
condition effect on the stability limit of cylindrical liquid bridges. 


3.2.1 Cylinder 

Johns and Narayanan determined the static stability limit of cylindrical liq¬ 
uid bridges with a pinned and a free contact line using perturbation techniques. 
Here, we derive the stability criteria using Myshkis’s variational method. Con¬ 
sider the cylindrical liquid bridge between two plates shown in Fig. H Plateau [l[ 
theoretically obtained the stability region A < 27r for cylindrical liquid bridges 
pinned at two equal coaxial discs (Fig. |4] (left)). In this section, the correspond¬ 
ing stability limit is obtained for cylindrical liquid bridges that are pinned to a 
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Figure 3: Equilibrium solution and existence region (above dash-dotted line) 
with respect to the canonical parameters; (a) constant-slenderness isocontours 
(label denote A), (b) constant-cylindrical volume isocontours (labels denote V), 
and (c) constant-slenderness and cylindrical volume isocontours. 
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Figure 4: Cylindrical liquid bridge: equilibrium surface is pinned to both discs 
(left); equilibrium surface is pinned at the upper disc and free to move on the 
lower plate (right). 


disc and free to move on a plate (Fig. |4] (right)). All the lengths are scaled with 
q as 


p=\q\r, ^ = qz, t = 


(31) 


with 



(32) 


used in Eqs. (fT^ - (ITCl) instead of . The solutions of Eqs. (fT^ - (ITCl) are 


wi(t) = sinr, W2(t) = cost, W3{t) = — 1 (33) 


for axisymmetric and 

mir) = T, W 5 {t) = 1 
for non-axisymmetric perturbations. These furnish 


(34) 


I?°(At) = —Ar sin Ar + 2(1 — cos At), (35) 


X°(At) 


D^{At) = At, 

At cos At — sin At 
—At sin At + 2(1 — cos At) ’ 

1 


x'(At) = - 


At ’ 


(36) 

(37) 

(38) 


where At = ti — tq and x = x/kl- The MSR boundary is determined 
by 11° (At) = 0 and D^{At) = 0 with respect to axisymmetric and non- 
axisymmetric perturbations, respectively. The first root of Eq. (1351) occurs 
where Ar^g^ = 27r, whereas Eq. (I36p has no non-trivial root. This implies 
that all cylindrical bridges with At > 27r are unstable to axisymmetric pertur¬ 
bations, irrespective of the contact-line condition at si. Note that the MSR of 
Fig. U (right) is equivalent to the stability region of Fig. U (left), and, therefore, 
the foregoing condition coincides with Plateau’s stability criterion. The first 
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non-trivial root oi = x and X^ = X inside the respective MSR identifies the 
stability region with respect to axisymmetric and non-axisymmetric perturba¬ 
tions, respectively. The former gives ~ 4.4934, whereas the latter has no 
non-trivial root. For cylindrical liquid bridges, A and At are equal; thus, the 
MSR and stability region can be summarized as A < 27r and A < 4.4934 where 
axisymmetric perturbations are the most dangerous. Note that the stability 
region of cylindrical liquid bridges with two free contact lines is A < tt l3£ 
which clearly indicates the destabilizing effect of free contact lines. 


3.2.2 Catenoid 

Here, we apply the same procedure as for cylindrical liquid bridges. Solving 
Eqs. (fT^ - (fTBll using the integral curve of Eq. (1^ gives 


wi(s) = 


\/P + I' 


W2(s) = 1 - 


+ 1 


ln(s -I- V+ 1), 


W 3 {s) = -- 


(39) 

(40) 

(41) 


W4{s) = - -I- 


2 V§ 


; ln(s -I- vi^ + l). 


W5{s) = 


+ 1 


(42) 


(43) 


Substituting Eqs. (IM]) - (H51) into Eqs. (fT71) . (fTSl) . (1^ . and (EH furnishes 
D^, and x^ as functions of sq and Si. These, unlike for cylinders, cannot 
generally be represented as functions of only As, implying that the stability 
of catenoids demands two independent parameters to be completely specified. 
This is consistent with the equilibrium solution discussed in section [3.11 

Erie et al. 2l| showed that catenoids pinned to two equal coaxial discs are 
unstable to axisymmetric perturbations when As/2 > 4.6395. We will deter¬ 
mine the stability region for catenoids with a free contact line as shown in Fig.[T] 
and demonstrate that they lose stability to axisymmetric perturbations. This 
can be accomplished by showing that = 0 and X^ = X have no non-trivial 
root (proved in Appendix . Accordingly, the MSR and stability region with 
respect to non-axisymmetric perturbations coincide with the existence region. 
Moreover, one can prove that = 0 has a non-trivial root only when 0c > 7r/2 
f Appendix Ell. This is clearly illustrated in Fig. [5l Here, D'^{6c,0d) is plotted 
in Fig. [5] (left) for 9c < 7r/2 and 9d € (tt — 0c,7r]. For a given 9c, —>■ 0“ 

as 9d ^ TT — 0c and —>■ — oo as —>■ tt. Thus, no sq can be found along 

the integral curve where vanishes, and the MSR spans the entire existence 
region. In contrast, for a given 9c > 7r/2, there exists a 9d (or sq) at which 
vanishes, as indicated in Fig. [5] (right). Here, ^ 0~ as 9d ^ tt — 9c and 
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Od Qd 

Figure 5: The effect of contact angle Oc on left panel: 9c = 10, 20, 30,50, 70° 
(blue, right to left) and 9c = 90° (red); right panel: 9c = 130,140,150,168,170° 
(blue, right to left) and 9c = 90° (red). 




Figure 6: Canonical phase diagram; (a) axisymmetric perturbations: constant- 
isocontours (thin-solid lines, labels denote x°), vanishingly small catenoids as 
—oo (thick-dash-dotted line), the MSR boundary as oo (thick-solid 
line), and the stability region boundary (dashed line); (b) non-axisymmetric 
perturbations: constant-x^ isocontours (thin-solid lines, labels denote x^)j and 
vanishingly small catenoids as x^ —oo (thick-dash-dotted line). 
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oo as TT. Note that 6 c — 168.75° is a special case because D° and 
dD^/dOd vanish simultaneously at dd — 162.07°. 

Figure IHl^a) shows the canonical phase diagram representing the MSR and 
stability region for axisyinmetric perturbations. The regions confined by Z'Z"mT'TU'Z' 
and Z'Z"mTZ' represent the MSR and stability region, respectively. The MSR 
boundary mT'TU' and stability-region boundary rriTZ' are determined, re¬ 
spectively, by Il°(so,si) = 0 and X°(so,si) = x(si). The meridian curve for 
catenoids corresponding to points on the MSR boundary satisfies ly = 0, and 
that corresponding to points on the stability region boundary satisfies A = 0. 

All the catenoids corresponding to points outside the MSR mUU'TT'm are, 
regardless of the contact-line condition at si, unstable to axisymmetric per¬ 
turbations. 11° and dD'^/ddc vanish simultaneously at T where {dc,6d) — 
(162.07°, 168.75°). Similarly, Zl° and dD^/dOd vanish simultaneously at T' 
where {6c, 6d) — (168.75°, 162.07°). Note that the MSR here is equivalent to 
the stability region of catenoids pinned to two unequal coaxial discs. Hence, 

Fig.m; a) also allows a comparison between two stability problems: (1) Catenoids 
pinned to a disc and free to move on a plate (Fig.[T]), and (2) catenoids pinned 
to two unequal coaxial discs with exactly the same sq and si. The region 
confined by mT'TU'Z'Tm represents the catenoids that are unstable to ax¬ 
isymmetric perturbations in the first problem, but stable in the second. lso-x° 
curves are thin solid black lines approaching Z'Z" {mT'TU') as ~oo 

(x° —>■ oo). Therefore, catenoids corresponding to points in the close vicin¬ 
ity of Z'Z" {mT'TU') are highly stable (unstable) since Xo X (Xo ^ x)- 
Figure [61(b) shows the canonical phase diagram representing the MSR and sta¬ 
bility region for non-axisymmetric perturbations. Here, iA^(sO)Si) = 0 and 
X^(sojSi) = x(si) have no non-trivial solution. Thus, the MSR and stability 
region coincide with the existence region, implying that catenoids are always 
stable with respect to non-axisymmetric perturbations. Iso-x^ curves are plot¬ 
ted as thin solid black lines approaching Z'Z" as x^ —oo. Note that x^ 

does not approach infinity for isocontours near the existence-region boundary 
Z'UZ". 

Figure ini also illustrates how the catenoid geometrical symmetry is reflected 
in its phase diagram. Catenoids that are pinned to two equal coaxial discs 
PH have equatorial symmetry, resulting in a one-dimensional phase digram in 
As. Even though catenoids bridging two unequal coaxial discs generally have 
no equatorial symmetry, and they require a two-dimensional phase diagram, 
a symmetric stability region can be constructed by choosing a proper set of 
parameters. For instance, one may choose the ratio of the lower and upper 
disc diameters K to represent the phase diagram (the second parameter can 
arbitrarily be selected). These catenoids are reflectively symmetric with respect 
to K. Clearly, inverting this ratio has no effect on the stability limit. Hence, the 
stability-region boundary must be invariant with respect to the transformation 
K = \/K. Alternatively, one can choose the dihedral angle that the catenoid 
forms with the upper disc 6 d and the lower one 6 c. The foregoing transformation 
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can equivalently be written 



which is why the MSR boundary in Fig. [S] (left) is symmetric with respect to 
the phase diagram minor diagonal described by Od = Oc- This is formally proved 
in Appendix [BJ Note that the stability-region boundary has no such symmetry 
since the contact line condition at i (see Fig.[T]) completely breaks the equatorial 
and reflective symmetries. 

Figure EKa) shows the existence-region boundary UTm in the favourable 
phase diagram (Om in Fig. [2]) and stability-region boundary Z'Tm. The curve 
UTm is the locus of points at which an iso-9d curve is tangent to an iso-0c curve. 
The curve corresponding to the existence-region boundary with respect to the 
favourable parameters UTm intersects the stability region boundary at T where 
{0c, 9d) — (162.07°, 168.75°). This is the point at which the slope of the MSR 
boundary is zero, as discussed for Fig. UKa). Figure EKa) clearly demonstrates 
that all the points on the segments UT and Tm correspond to unstable and 
stable catenoids, respectively. Figure EKb) shows isocontours in the canonical 
phase diagrams. Selecting two variables among A, V, 9c, and 9d, this figure 
describes the equilibrium solution and stability of the corresponding catenoids. 
Consider the point (A, V) ~ (1.0518,1.0264) in the favourable diagram, for ex¬ 
ample. It lies at the intersection of 9c = 130° and 9d = 110°. This can be 
located in the canonical phase diagram, as shown in Fig.EJb). The correspond¬ 
ing point A is inside the stability region, indicating that the respective catenoid 
is stable. Furthermore, the second equilibrium solution can be determined by 
identifying the other intersection point of the same iso-A and iso-R. This oc¬ 
curs at A', where {9c, 9d) — (177.08°, 170.1°). The second solution lies outside 
the stability region, which corresponds to an unstable catenoid. Consider the 
point (A, V) ~ (0.5555,0.1742) in the favourable diagram. The two equilibrium 
solutions are represented in Fig. EKb) by D at {9c, 9d) = (130°, 170°) and D' at 
{9c, 9d) — (175.32°, 174.96°). Both equilibrium solutions lie outside the stability 
region and correspond to unstable catenoids. 

Volume-controlled catenoids (constrained) are exposed to a smaller set of 
disturbances than pressure-controlled catenoids (unconstrained). Therefore, the 
former is more constrained and stable. Using Maddocks’ theorems 0 , Akbari 
et al. [HI showed that catenoids lose stability with respect to constant-pressure 
perturbations at turning points in A along iso-0c curves (see Fig.EKb)). Accord¬ 
ingly, the upper and lower segments of these isocontours respectively correspond 
to stable and unstable catenoids. Figure EKb) compares constant-volume and 
constant-pressure stability regions. Here, the region confined between the long- 
dashed and dashed-dotted lines is the constant-pressure stability region and 
that between the thick-short-dashed and dashed-dotted lines is the constant- 
volume stability region. As expected, the constant-pressure stability region is 
completely contained inside the constant-volume stability region, indicating that 
constrained catenoids are more stable than unconstrained catenoids. 
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Figure 7: Existence region in the canonical phase diagram, existence-region 
boundary (dash-dotted for canonical and thin-short-dashed for favourable), 
stability-region boundary with respect to volume-controlled (thick-short- 
dashed) and pressure-controlled (long-dashed) perturbations; (a) stability re¬ 
gion in the canonical phase diagram with (b) the respective isocontours (the 
same as Fig. ED- 
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4 Concluding remarks 


We have examined the equilibrium and stability of catenoids bridging a circular 
disc and plate where the equilibrium surface is pinned at one contact line to the 
disc edge with the other free to move on the plate. Drawing on the second vari¬ 
ation of potential energy, the existence, maximal stability, and stability regions 
were analytically determined. The equilibrium solution multiplicity subject to 
various constraints was discussed in detail. The results showed that all catenoids 
are stable with respect to non-axisymmetric perturbations; for a fixed contact 
angle, there exists a critical volume below which catenoids are unstable to ax- 
isymmetric perturbations. The canonical phase diagram furnishes a one-to-one 
correspondence between points in the existence region and equilibrium solutions 
where the stability-region boundary separates the points corresponding to stable 
catenoids from those corresponding to unstable ones. No such correspondence 
can be established in the favourable phase diagram. Furthermore, the canonical 
phase diagram conveniently demonstrates how the catenoid geometrical sym¬ 
metry affects the stability regions. For example, the maximal stability region 
symmetry with respect to the phase diagram minor diagonal indicates the re¬ 
flective symmetry (with respect to the ratio of lower and upper disc diameters) 
of catenoids with two pinned contact lines. Moreover, the asymmetric shape 
of the stability region shows how a catenoid free contact line with a substrate 
breaks the equatorial and reflective symmetries. The stability limit presented 
here is a limiting case for the minimum volume stability limit of liquid bridges 
when the mean curvature approaches zero The static stability limits are 


useful for predicting the transition of the time scale from the quasi-static to the 
intermediate phases of contact-drop dispensing. 
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Appendix A Stability with respect to non-axisymmetric 
perturbations 

We shall prove that all catenoids are stable with respect to non-axisymmetric 
perturbations. First, we show that D^(so,Si) = 0 has no non-trivial root. 

Substituting Eqs. (I5^ - (H5]) into Eq. (1^ yields 

si\Js\ 1 + ln(si -|- + 1) = So1 + ln(so + -I- 1), (A.l) 

where we seek an sq £ (—oo,si) for a fixed si. Equation (lA.ll) is obviously 
satisfied for the trivial solution sq = si. Denoting the left-hand side by /(si). 
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Eq. (lA.ll) can be rewritten as 



Sq + 1) —/(si) — 0 (A.2) 


with 



(A.3) 


All the terms on the right-hand side of Eq. (IA.3|) are positive, indicating that 


dg/dso > 0 on (—oo,si). The continuity of ^(so) implies that 5 is a mono- 
tonically increasing function such that g —> —00 as sq —> —00 and g —>■ 0 “ as 
So —si- Therefore, sq = si is the only solution of Eq. (IA.2I) . 

Next, we prove that x^(so,si) — x(si) = 0 has no non-trivial root. Substi¬ 
tuting Eqs. (I5^ - (H51) into Eq. (IT^ results in 



Note that the denominator of the fraction in the square bracket is non-zero for 

50 € (—oo,si). Two casesmust be considered separately: (1) Si 0 and (2) 

51 —>■ 0. The hrst leads to 



(A.5) 


A pair (so,Si) that satisfies Eq. (IA.4I) must also satisfy Eq. (IA.5I) . Equa¬ 
tion (IA.5I) clearly shows that Eq. (IA.4I) has no non-trivial root since it is inde¬ 
pendent of So. In addition, no si can be found that satisfies Eq. (IA.5|) because 
the left-hand side is always greater than zero. Erom the second case, 


(A. 6 ) 


So \/so + 1 ~ ^ So 


which holds only when so —>■ — 00 . This indicates that, for any given contact 
angle, only infinitely large catenoids may lose stability to non-axisymmetric 
perturbations, and the stability-region boundary coincides with the existence- 
region boundary in the canonical phase diagram. 

One can apply the same procedure as the two previous cases to prove that 
I?°(so,si) = 0 has no non-trivial root for fy < Tr/2. However, the expressions 
are cumbersome and the analysis is tedious. We only demonstrate the limiting 
behaviour discussed in section 15.2.21 The Taylor-series expansion of is used 
for the small-interface limit: 



(A.7) 


where e = si — Sq- It follows that 


lim D° = 0-. 


(A. 8 ) 


e^0+ 
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In the other limit, where catenoids are infinitely large, one can show that 

lim = sign(si) x oo, (A.9) 

Sq—^ —OO 

implying that there is at least one sq G (— 00 , si) at which D°(so,si) = 0 for 
Oc > 7 r/ 2 . These limits are clearly illustrated in Fig. [S] 


Appendix B Symmetry of Z)-functions 


We prove that and are symmetric with respect to the canonical phase 
diagram minor diagonal. Consider the following transformation 


r So = -Si 
\ Si = -So 


(B.l) 


which is equivalent to Eq. dg. Given that rci, wa are odd and r, W 2 , W 3 , W 5 
are even functions, we have 


£'°(so,Si) = 


-wi(si) W 2 {h) 

-wi(so) W 2 (so) 

— fwids fw2ds 


wsisi) 
Wsiso) 
fs7 ^^3ds 


D^(so,Si) 


-W4(si) W5(si) 

-W4(so) W5(so) 


(B.2) 


(B.3) 


Taking the determinant row exchange rules into consideration, it follows that 


A’°(so,si) = i 4 °(so,Si), 

(B.4) 

£•^(50,51) = i 4 ^(so,Si). 

(B.5) 


This completes the proof. 
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